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Macromolecules occupy between 34 and 44% of the cell cytoplasm, about half the maximum pack- 
ing density of spheres in three dimension. Yet, there is no clear understanding of what is special 
about this value. To address this fundamental question we investigate the effect of macromolecular 
crowding on cell metabolism. We develop a cell scale flux balance model capturing the main features 
of cell metabolism at different nutrient uptakes and macromolecular densities. Using this model we 
show there are two metabolic regimes at low and high nutrient uptakes. The latter regime is charac- 
terized by an optimal cytoplasmatic density where the increase of reaction rates by confinement and 
the decrease by diffusion slow-down balance. More important, the predicted optimal density is in 
the range of the experimentally determined density of E. coU. We conclude that cells have evolved 
to a cytoplasmatic density resulting in the maximum metabolic rate given the nutrient availability 
and macromolecular crowding effects and report a flux balance model accounting for its effect. 



Macromolecular crowding affects the rate of biochem- 
ical reactions [l]]. It tends to increase reaction rates by 
increasing enzyme concentrations and to decrease reac- 
tion rates by reducing the diffusion coefRcient of metabo- 
lites. The competition between these two factors results 
in a maximum reaction rate at intermediate crowding 
agent concentrations [l[. Yet, it remains to be addressed 
whether this observation is true at the level of the whole 
cell metabolism. The macromolecular volume fraction 
of the E. coli cytoplasm is in the range 0.34-0.44 0, 
half the maximum packing density of spheres in three di- 
mension. This observation indicates cell metabolism is 
operating under the crowding conditions created by its 
macromolecular components. 

To investigate the effect of macromolecular crowding 
in the overall cell metabolism we focus on the model 
schematically represented in Fig. [TJ (i) From the 
perspective of cell metabolism alone the cytoplasm is 
composed of metabolites and enzymes, the latter in- 
cluding metabolic enzymes, ribosomes and any other 
macromolecule catalyzing metabolic/biosynthetic pro- 
cesses, (ii) Metabolites are relatively small compared to 
enzymes and we neglect their contribution to crowding, 
(iii) The active site of most metabolic enzymes is rel- 
atively small compared to the whole enzyme. For all 
practical purposes the inert enzyme region is equivalent 
to an inert crowding agent. The inert region reduces the 
volume available to all solutes and, from a point of view 
of diffusion, collisions between metabolites and the inert 
region of enzymes are equivalent to collisions between 
the metabolites and inert crowding agents. Therefore, 
we model the contribution of active sites and the in- 
ert mass independently. We assume that each enzyme 
molecule contributes as two different fictitious quasi- 
molecules: one representing the enzyme active site (or 
union of active sites) and the other the inert enzyme re- 
gion (Fig. [Th). The active site quasi-molecule is assumed 



to occupy a relatively small volume and its contribution 
to crowding is neglected. Its diffusion coefficient is as- 
sumed to be, however, the same as that for the original 
enzyme. The quasi-molecule representing the inert en- 
zyme region is modeled as an inert crowding agent, with 
a size and diffusion coefficient equal to that of the corre- 
sponding enzyme. Because of their relatively larger sizes, 
enzymes have smaller diffusion coefficients than metabo- 
lites. We thus approximate the relative diffusion coeffi- 
cient between enzymes and metabolites by the diffusion 
coefficient of the metabolites. 

(iv) In a crowded media concentrations are effectively 
higher because part of the volume is occupied by the 
crowding agents. The ratio between the effective concen- 
tration C and the concentration in an ideal solution Co 
is denoted by the activity coefficient 7c = C/Cq- When 
the solute interacts with the crowding agents exclusively 
via steric repulsion = V/Va 0, where V is the total 
volume and Va is the volume available to the solute. The 
volume available to a solute can be approximated by the 
total volume minus the volume occupied by all crowding 
agents, resulting in 



where v is the macromolecular volume fraction. 

(v) The diffusion coefficient of a trace particle in a 
crowding media is also affected by the concentration of 
crowding agents. In general the diffusion coefficient is 
given by D = JdDq, where Dq is the diffusion coefficient 
in aqueous solution and 7d is a correction factor. To 
quantify the impact of crowding on metabolites diffusion 
we use the empirical exponential law 

7D = e-"'' , (2) 
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FIG. 1: Metabolic model of the cytoplasm, a) Schematic 
representation of the cytoplasm from the perspective of cell 
metabolism. The solid circles represent small metabolites. 
The dashed circles with attached small empty circles repre- 
sent metabolic enzymes, the former representing their inert 
region and the latter the active site, b) We model an enzyme 
of molecular volume v and diffusion coefficient D as two dif- 
ferent quasi-molecules. One representing the active site, with 
zero (or negligible) molecular volume and the same diffusion 
coefficient. The other representing the enzyme inert region, 
with same molecular volume and diffusion coefficient as the 
whole enzyme. 

where the exponent a is an empirical parameter 0, [^, 0, 
0|. The closest estimate for the cytoplasm comes from 
the experimental report a = 5.8 for fibroblast cells fHj. 

(vi) The kinetic models describing the rate of biochem- 
ical reactions as a function of the concentration of react- 
ing metabolites can be quite complex and are unknown 
for most metabolic reactions. For example, consider the 
simple case where a substrate S is irreversibly trans- 
formed into the product P catalyzed by the enzyme E, 
with concentrations 5*, P, and respectively. The rate 
of this reaction is given by the Michaelis-Menten model 
R — k2SE /{Km + S), where ^2 is the rate of the conver- 
sion of the intermediate enzyme-substrate complex ES 
into the product ES-^P, Km — 4- fc2)/fci is the 

Michaelis-Menten or half-saturation constant, and 
and ki are the rate of the intermediate steps ES^E-|-S 
and E+S^ES, respectively. The step E-hS^ES is dif- 
fusion limited and fci — AnDa, where D is the sub- 
strate diffusion coefficient and a is the effective size of 
the enzyme active site. Diffusion limited regime: When 
S <C Km the reaction rate can be approximated by 
R Ki [fc2/(fc-i -t- k2)]kiSE (X DE and therefore the over- 



all reaction is diffusion limited. In this limit most active 
sites are free and the reaction rate is limited by the rate 
of encounter of the enzyme active site and the substrate. 
Saturation regime: When S ^ Km the reaction rate is 
approximated by i? w k2E. In this case most active sites 
are occupied and the reaction rate is limited by the chem- 
ical step ES^P. The situation becomes more complicated 
for reactions involving more than one substrate, because 
of reversibility and other factors. Nevertheless, in gen- 
eral these two limiting scenarios, diffusion limited and 
saturation persist. As a first approximation we therefore 
assume that biochemical reactions are divided into two 
groups: a set C of diffusion limited reactions and a set S 
of reactions at saturation, with reaction rates given by 

_ J giDid , for i e £ , , 

~ \ k,a for i e 5 ^ ' 

where gt is a model parameter containing all other contri- 
butions in the diffusion limited regime (e.g. metabolite 
concentrations) and ki is the rate of reaction i in the 
saturation limit. 

(vi) At the level of cell metabolism we use a steady 
state or flux balance model In the steady state the 
consumption and production of each metabolite balance 

^S,,R, = 0, (4) 

i 

where i = 1, . . . ,n as above is an index over reactions, 
J = 1, . . . , TO is an index over the metabolites, and Sji 
is the the stoichiometric coefficient of metabolite j in 
reaction i Q . To account for the potential existence of a 
limiting nutrient we label by i = 1 a reaction representing 
the nutrient uptake and assume 

Ri<U , (5) 

where U denotes the maximum uptake of the limiting 
nutrient. We also label hy i — n the metabolic objective 
or biomass vector, an effective reaction with a nonzero 
stoichiometric coefficient for each metabolite the cell pro- 
duces and magnitude given by its relative ratio. Thus i?„ 
is our measure of metabolic rate. 

(vii) Given a flux distribution Ri we now calculate the 
volume fraction occupied by enzymes. The metabolic en- 
zymes are characterized by their concentration Coi, oc- 
cupied volume fraction Vi, molar mass fii, specific volume 
Vi and kinetic model ([3]). The occupied volume fraction 
is related to the enzyme concentration through the equa- 
tion Vi = fj^iViCoi. In turn the enzyme concentration is 
related to the reaction rate through Putting these 
two relationships together, recalling that Ci = jcCoi and 
Di = joDoi, and using ^ and we obtain 
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FIG. 2: Calculation of the optimal volume fraction. 

The left (solid line) and right (dashed lines) hand side of ((T]) 
for different values of U, assuming S/L = 2 and UL = 0.03, 
0.057 and 0.08 from bottom to top. The intersection between 
these lines determines the optimal volume fraction v* . For 
the upper dashed line, representing a larger U , there is no 
solution. 



FIG. 3: Optimal volume fraction. The optimal volume 
fraction for a = 5.8 as a function of the ratio S/L (solid 
line). The shadowed area was obtained using as input the 
volume fraction range 0.34-0.44 for E. coli 0] and computing 
the associated S/L range from (O. The dashed line represents 
our prediction v* — 0.37 assuming S/L = 2>. 
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At this point we add up the volume fraction occupied by 
enzymes v = u^, which after some algebra results in 



l-v 



U[L{r) + S{ry 



L(r) 



S{r) 



(7) 



(8) 



where — Ri/U are the reaction rates in units of U . 

Under the model (i)-(vii), the optimal reaction rates 
i — l,...,n, and volume fraction v are obtained 
maximizing the metabolic objective subject to the 
flux balance (|4]), uptake capacity (O and solvent capac- 
ity ([7])- ([5]) constraint. The solution of this optimization 
problem results in the following regimes. 

Nutrient limited: There is a threshold uptake capacity 
Uc such that for U < Uc the optimal reaction rates, de- 
noted by r^^\ are those maximizing the metabolic objec- 
tive Rn subject to the flux balance ^ and uptake capac- 
ity ([5]) constraint, r^^"^ are independent of U, and there 
is a volume fraction v satisfying ([7]) with L — L{A^^) 
and S ~ S{A^'>). The existence of this threshold uptake 
capacity is derived from the analysis of solutions to ([7]) 
with respect to v, given L and S (Fig. [5]). The left hand 
side has a maximum at 



On the other hand the right hand side is a decreasing 
function of u, starting from U{L + S) at w = and end- 
ing at U{L + Se^^"^) at v — Vm, where Vm denotes the 
maximum packing density of the cytoplasmatic compo- 
nents. For small values of U the two curves intercept and 
([7| has a solution v = v{U L, S/L), which is an increasing 
function of U. There is no solution, however, for large 
values of U . The maximum values of U and v where a 
solution exists, denoted by Uc and v* , are obtained when 
the tangents of the left and right hand side of ([7]) are 
equal as well. The simultaneous solution of these two 
equations results in an explicit relation between Uc and 
V* , and a transcendental equation for v* parametrized by 
a and S/L. For the intracellular value a ~ 5.8 we com- 
pute the maximum cell density as a function of S/L (Fig. 
[3]). It increases monotonically starting from vq w 0.22 ^ 
when all reactions are diffusion limited. 

Space limited: For U > Uc 'Vfe cannot further increase 
the volume fraction beyond v* and ([7]) becomes a con- 
straint on the reaction rates. In this regime the optimal 
reaction rates are obtained maximizing the metabolic ob- 
jective Rn subject to the flux balance dU, the uptake 
capacity (O, and the solvent capacity constraint 



atRi = 1 



giDoi 



for i G L 
for i G S* 



(10) 



(11) 




< a < 4 
a > 4 . 



An equation similar to (|10|1 has been introduced before 
^g^ under the name of macromolecular crowding or solvent 
capacity constraint The coefflcients have been 



4 



named crowding coefficients as they quantify how much 
the reactions contribute to the crowding of the cyto- 
plasm. Equation pT|) now introduces corrections to the 
previous calculations , making a more precise account- 
ing for the effect of macromolecular crowding on reactions 
rates. 

To provide a quantitative prediction for the optimal 
volume fraction we need an estimate of SjL. Given that 
the kinetic parameters involved in these calculations are 
unknown for most biochemical reactions we cannot make 
a precise calculation at this point. Recent metabolite 
concentration measurements for E. coli indicate that 
83% of the reactions have S > Km, i-e. there are about 
three times more reactions in the saturation regime. As- 
suming that the values of gtDoi and ki are not signifi- 
cantly different we would conclude that S is about three 
times larger than L. In this case our calculations pre- 
dict an optimal volume fraction of 0.37 (see Fig. [3]). 
Experimental estimates of the macromolecular volume 
fraction of the E. coli cytoplasm indicate the lower and 
upper bounds 0.34 and 0.44 respectively 0, containing 
the crude estimate 0.37 predicted above. A more conser- 
vative prediction is the interval between the lower bound 
0.22 when all reactions are diffusion limited and the up- 
per bound given by the maximum packing density Vm 
when all reactions are at saturation. The latter can reach 
values as high as 0.8 when the particles have variable size 
llj. This conservative range also contains the experi- 
mental report between 0.34 and 0.44 for E. coli 

In summary, there are two different metabolic regimes, 
nutrient limited and nutrient limited with a solvent ca- 
pacity constraint. In the nutrient limited regime the 
metabolic rate is constrained by the stoichiometry of the 
set of biochemical reactions and the maximum uptake 
rate of the limited nutrient. In the nutrient limited with 
a solvent capacity constraint the metabolic rate is in ad- 
dition constrained by macromolecular crowding effects. 
In particular there is an optimal macromolecular volume 
fraction resulting in the maximum metabolic objective. 
The optimal volume fraction is a function of the diffu- 
sion related exponent a and the ratio of reactions in a 
diffusion limited and saturation regimes. It interpolates 
between v* ~ 0.22 when all reactions are diffusion lim- 
ited and the maximum packing density when all reactions 
are at saturation. Finally, we can write a general flux 
balance model that applies for both regimes: maximize 
the metabolic objective i?„ subject to the flux balance 
J2i SjiRi = 0, uptake capacity Ri < U and solvent ca- 
pacity UiRi < 1 constraint. The latter inequality is 
satisfied with the less than sign for U < Uc and with the 
equal sign for U > Uc- 

We have made several assumptions that could affect 
the obtained results. The equations ^ and ^ charac- 
terizing the impact of macromolecular crowding on ef- 
fective enzyme concentrations and metabolite diffusion 



coefficients may require further corrections when the oc- 
cupied volume fraction gets closer to the maximum pack- 
ing density. In the intracellular media the typical occu- 
pied volume fraction is about half the maximum packing 
density 0.74 of spheres in three dimensions. Therefore, 
we expect equations ([T]) and ^ to be sufficiently good 
approximations. The kinetic models cannot always be 
approximated by the extreme cases of diffusion limited 
and saturation regimes. As a consequence the plot in 
Fig. [3] may be slightly different. 

Taking together our results indicate that at high 
metabolic rates there is an optimal intracellular density 
where the increase of reaction rates by confinement and 
the decrease by diffusion slow-down balance. Since it is 
the optimal density resulting in the maximum metabolic 
rate, an increase of the density of enzymes beyond this 
optimal value will result in a decrease of the metabolic 
rate. Although this may sound counter intuitive, it fol- 
lows from the fact that beyond the optimal density the 
slow-down of diffusion starts to dominate, diminishing 
the overall metabolic rate. More important, the exper- 
imentally determined density of E. coli is in the range 
predicted from our model. We thus claim that cells have 
evolved an intracellular density that results in maximum 
metabolic capabilities given the macromolecular crowd- 
ing effects. 
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